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Abstract. In this work, we develop a framework based on piece wizc 
Bezier curves to plane shapes deformation and we apply it to shape opti- 
mization problems. We describe a general setting and some general result 
to reduce the study of a shape optimization problem to a finite dimen- 
sional problem of integration of a special type of vector field. We show a 
practical problem where this approach leads to efficient algorithms. 

In all the text below, £^ = IR^. In this text, we will define a set of manifolds, 
each point of such a manifold is a parametrized curves in E. 

1 Bezier curves 

Bezier curves are usual objects in Computer Aided Geometric Design (CAGD) 
and have natural and straightforward generalization for surfaces and higher di- 
mension geometrical objects. We focus here on curves even if a lot of results have 
natural generalization in higher dimension. This section has aim to fix notation 
and make the paper as self contained as possible. 

1.1 Beisic definitions 

Given Po,Pi,. . . , Pn e E, we define: 

B{{Po,...,PD),t)={l-t)B{{Po,...,PD-l),t)+tB{{P,,...,PD),t) 

with B ((P) ,t) = P for every P G E. The associated Bezier curve is 

{B {{Po, Pd) ,t) \ t G 0,1} and the list {Po,...,Pd) is called the control 

polygon and the points Pq, . . . , Pd arc called the control points. 

This process associates to every set of points a parametrized curve. It is a 
polynomial parametrized curve and its degree is bounded : 

Proposition 1. Let Pq, ■ . ■ , Pd € E, then B {{Po, ■ ■ ■ , Pd ) ,t) is a polynomial 
parametrization and its coordinates have degree at most D. 



1.2 Bernstein's polynomials 

Definition 1. Let D be an integer and i G {0, . . . ,D}, we define the Bernstein 
polynomial bi £, (t) :— (^) (1 — i) * t^ . 

Notation 1 We denote IEl[i]jj the set of polynomial of degree less or equal to 
D. The set R [t]jj has a natural '^-vector space structure, its dimension is Z? + 1 
and is a basis of this vector space. 

Proposition 2. The set {bo^jj, . . . jbjj^D} is a basis o/R[t]^. 

N 

Proposition 3. Let Pq, . . . , e E, then B {{Pq, . . . , Pd) ,t) = J^Ptbi.oit) 

1=0 

for all t e [0, 1]. 

Corollary 1. Every polynomially parametrized curve can be represented as a 
Bezier curve. 



1.3 Interpolation 

Since a Bezier curve of degree D is defined using D + 1 control points, one can 
hope to associate D + \ control points from a sampling of + 1 points on a 
curve. The following result shows that this is possible. But in fact, we do not 
have one Bezier curve of degree D but many ones. Each such curve is associated 
to a particular sampling of the parameter interval [0, 1]. 

Proposition 4. Let Mq, . . . , Md e E, then there exists Bezier curves of degree 
D passing through these points. 

Lemma 1. Let to = < ti < ■ ■ ■ < t^ = 1, then there exists one and only one 
Bezier curve B {{Po,...,Pd) ,t) of degree D such that B ({Pq, . . . , Pd) ^t,) = 
Af„Vi e {0,...,!?}. 

Proof. Denote M the 2 x (D + 1) matrix built with the coordinate of Mi as i^^ 
row, i.e. M = (Mq, . . . , M^Y, and denote P the 2 x {D + 1) matrix built with 
the coordinate of Pi as i*"^ row, i.e. P = {Pq, . . . , PdY ■ We consider the following 
matrix associated to t — (^to, . . . ,t]j^ : 



Btn — 



( bo,D (0) (0) • • • bD,D (0) \ 

bQ,D (tl) bi^D (ti) ■ ■ ■ bD,D {ti 



(1) 



V 60,13(1) Knil) ■■■ bD,D{l) J 

The matrix of equation [l] is invertible (it is the Vandermonde matrix express in 
the Bernstein basis) and clearly if P is such that Bt^nP = A/, then B {[Pq, . . . , Pd] 
give the wanted curve for the proof of the lemma. 

Remark that once t is know one can compute B^^ once for all and that is 
possible to take advantage of its Vandermonde-like structure in order to improve 
the cost of the multiplication of a vector by Bt^D- Generally, we use a regular 
subdivision {ti = j^) but there are more suitable choices in regard of the stability 
of the computation. 



2 Piecewize Bezier curves 



2.1 Basics on piecewize Bezier curves 

Let Po,Oi • • • ) Po,N and Pifl, ■ ■ ■ , Pi,d S E such that Pq^n = Pi,o, we define the 
following parametrization of a curve: 

f [0,1]^^; 

,^ fB((Po,o,...,Po,Z3),2i)fortG [0,1/2] 
\ B((Fi,o,...,Fi,i)),2i-l)forte [1/2,1] 

When N = D we say that this parametrization is uniform with respect to the 
degree and often, we simply say uniform when it does not introduce ambiguity. 
The curves parametrized by B (( Pq^Qi • • • j Po,d ) , i) and B (( Pi o, • • • , Pi,d ) , i) 
are called the patches of C = P ([0, 1]). The set of control points of the patches 
of C are called the control points of C. This is a continuous curve. 

More generally, if Po,o, • • • , Po,Ni , Pi,o, • ■ • , Pi,Ar2 , ■ ■ ■ , Pifl, ■ ■ ■ , Pl,Ni € E, such 
that Pi^Ni = Pi+1,0 for alH G {0, . . . , i — 1}, we define: 

r (((-Po,o, • • • , Po,d) , • • • , (Pat.o, • • • , Pn,d )),t)=B ({Pi,o, . . . , Pi,D,) ,j^ + {l + l)t 



for t G 



/ + 1 

(2) 

(1+1)7(1+1) and for all i G {0,...,/} . This defines a continuous 

parametrization. The curves parametrized by B {{Pi,o, ■ ■ ■ ,Pi,Di ) ,t) are called 
the patches of C = P ([0, 1]). Furthermore, if Pi^Di = Po,o we say that the curve 
C is closed or that it is a loop. 

We denote Bn,d the set of uniform piecewize Bezier curves built from A'^ 
patches of degree D. This clearly a finite dimensional subvariety of ([0, 1] , E) 
as the image of the following map: 



N,D 



-\{{Pij,j = O...D),i = O...N)^r{{{Pij,j = O...D),i = O...N),t) 

(3) 

Clearly, ^n,d is onto fi-om (P^^+i)^+^ to Bn,d C C° ([0, 1],E). It not very 
difiicult to check that ^n,d is almost always one-to-one from to 

Bn,d- So, tf'jv,D is almost everywhere a difFeomorphism between (^E^^^'j^'^^ 
and Bn,d- This embed B^^d with a manifold structure (even a submanifold 
structure in C°([0,1],E)). 

The density of polynomials in the set of continuous functions imply that for 
each <P : [0, 1] — > E there exists (P„ (t))„pw such that lim 1|^ - P„|L = 0, in 
a way that considering Bezier curves is not a drastic restriction. 

2.2 Sampling map and retraction to ^N,n 

Definition 2. Let to — < ti < . . . < to = 1, we denote t — {to, . . . , to) the 
associated subdivision of [0, 1], then we define the sampling map St ■ Bi d — ^ 

E'^+^ bySt{r) = {r{to),...,r{tD)). 



Proposition 5. The following diagram is commutative: 

ED+1 

and ^N.D is an invertible linear isomorphism between Bi^d — Im {^i,d) o.nd 
E^^^ and its inverse is S'^^j = o St- 

N D 

Proof. Let ^ (i) = E EPJ■^b^.D (t) i.e. F = ^n.d {{Po.q, • • • , ^^0,^) , ■ ■ • , [Pn.o, • ■ • , Pn.d)), 

j=li=0 

I \ 

then clearly St {F) = Bt.oP where P = : , and so St o ^Fn.d {P) = 

[Pkoj 

Bt,N (P)- The remainder of the theorem is a consequence of the fact that Bt^N 
is a linear isomorphism. 

Proposition 6. Let ti,o = < ii^i < ••• < ti,D = 1/^^ = ^2,0 < h.i < 
■ ■ ■ < t2.D — = ^3,0 < • • • < tiq^D = 1, we denote t = (ti, . . . , t^) where 
U — {to,i, ■ ■ ■ ,'tD,i) and we define the sampling map St.N ■ Bm,d — > [E^)^ by 
St,N (-T) = X • • • X 5t„ (F) = (r (ti^o) , . . . , r {tM,D) ) • Then St^N is a linear 
isomorphism between Bm.d o,nd (E^^^^ ^ . 

Proof. It is a simple consequence of the fact that a cartesian product of iso- 
mophisms is an isomorphisms. The inverse map is the cartesian product of the 
inverse of the component maps. 

The proposition [6] is important since it allows to give to Bm,d a vector space 
structure isomorphic to (and so, of finite dimension). For instance, 

it allows to transport distance and so on in Bpf jj. 

In fact, we focus here into a speciale type of sampling. We consider such an 
sampling where ti ^ = and ti^u = and tij = ti^ + We will call this 
kind of sampling a regular sampling and we will omit the subscript t everywhere 
using these samplings. We use these sampling to simply the presentation, but 
all the results has equivalent statements with general sampling. Representing 
each patch by its control polygon, the matrix of St,N is N times the cartesian 
production of the map with itself: Bi £, x • • • x Bi q. This gives us an easy 
way to solve the following interpolation problem. 



Problem 1. Given 7\fo,oi ■ • • i -^^o,Di • • ■ 7 -^^Ar,0; ■ • • : ^'^a^.d G B, find F e Bm,d 
such that Sn {F) - 



Proposition 7. The solution of problem^is given by the image by 'I'n.d of: 

(5) 



^i,D \ I ^"0,0 




The proposition m implies that Xt.D — B^d ° St '■ Bn.d — > E^'^^ is such 
that tf'i.D °Xt,D = lagD+i. It is easy to extend this resuh to ^n.d using Bfq d = 
Bi^D X • • • X Bi^D satisfying B^'^jj = B~^jj x • • • x B~^jj. 

This approach allows to project any element of C° ([0, 1] , E) on Bm.d using 
5jv. Let .1 G C ([0,1] ,£;), then denoting M = (r(0),r(^),.. .,r(^^),r(l))* 
we have P = SJ^^ {M) £ Bn,d is such that ^n,d (P) — B {P, t) coincides with 
yl([0, 1]) on at least (D + 1) points counted with multiplicities on each patch. 
This is only the fact that xt,D can be extend to C° ([0, 1] , i?)- 

The main claim is that instead of working directly with Bn.d, it is easier 
to work on the "set of control polygons", namely E^'^^ using sampling and 
interpolation giving linear isomorphism between control polygons and sampling 
points on the curves. In what follows, we will always take this point of view. 



2.3 Tangent space TBn,d and deformation of curve 

Recall that ^n,d define a linear isomorphism between the "space of control poly- 
gons 

(E^+Y+' and the space of piecewize Bezier curves Bn,d- We already 

saw that for any 7 (t) e Bnm then P G such that Wn,d (P) = 7 {t) 

is given by B^^j^ o Sm (7) . This give the following proposition: 

Proposition 8. We have that T^n.d ■ T (£'^+i)^^^ — > TBn.d is such that 
from any 7 G Bn,d we have T<P^^jj (7) : T^Bn,d — > '^xn,d{i) {E^^^)^^^ is 
given by T<I'n,d {xnm (7))"^ (e) B^^jj o Sn (e) = Xn,d (e) for any e (t) G 
TjBn.d ttnd this is a linear isomorphism. 

An element of e (i) G T^Bm.d is called a deformation curve. In fact, this 
proposition allows to express, given a piecewize Bezier curve and a deformation, 
how to deform its control polygon. This is an essential step proving that manip- 
ulating piecewize Bezier curve, it is enough to manipulate its control polygon. 
This is the object of the following lemma. 

Lemma 2. Let P G {E^+^)^^^ , j (t) = 'Pn,d{P) = B{P,t) G Bn.d and 
e{t) G T~^Bn,d, then: 



i. e (t) = IfN^D {XN,D (e))- 

a. j {t)+e {t) = ^N^D {P + XN,D {£))- 



This lemma explain how to lift a deformation from the space of curves to the 
space of control polygons. The vector space structure of both the space of control 
polygons 

^ED+i)"+' and of piecewize Bezier curves Bn,d allows to avoid the 
use of computationally difhcult concept as exponential map between manifold 
and its tangent space and so on. This structure has also to define a simple notion 
of distance between two such curves. 



3 Applications to shape optimization 

In this section, we show how the preceding formalism can be exploited in the con- 
text of shape optimization. An application to a problem of image segmentation 
is presented to illustrate our purpose. 



3.1 Shape optimization problem 

A shape optimization problem consists in, given a set of admissible shapes A 
and a functional F : A ^ find a shape a ^ A such that for all other 
shapes f3 G A, we have F (a) ^ F {(3). Generally, one try to give to the space of 
admissible shape a structure of manifold in a way to be able to compute a "shape 
gradient" , WF{/3), expressing the evolution of the criterium F with respect 
to a deformation of the shape (3. It is to say that VF {(3) associates to every 
point M G /3 a deformation vector VF (/3) (M) e TmE. The computation of 
such a gradient can require sophisticated computation since very often, even the 
computation of the criterium itself require to solve a system partial differential 
equations. Many problem can be expressed as a shape optimization problem. 
Classical approach to solve this kind of problem is to use VF(/3), when it is 
computable, in a gradient method to find a local minimum. 

To keep the presentation as simple as possible, we focus on geometric opti- 
misation, i.e. the topology of the shape is fixed, in the case where the frontier 
of the admissible shapes are continuous Jordan curves. But the framework pre- 
sented here can be extended to topological optimization as it is shown in [5^ for 
a special application on a problem of image segmentation. The case treated here 
received attention because of its deep links with images segmentation and shape 
recognition (see [5 3 2J for instance). 

We denote Cj ([0, 1] , F) the set of function parametrizing a Jordan curve and 
B% jj = {7 e Bn,d I 7 (t) = 7 (s) withs ^ t ^ {t = Q ands = 1) or {t = lands = 0)}. 
We have B% C C°j ([0, l],E). We denote: 



H 



N,D = {m.j,3 ^O...D),i^O...N)€ {E^'f I Po.o - Pn,d] 



which is a linear subspace of (£'^+^)^^^. We then denote 'I'f^ ^ = 'i'N,D \hn.d ■ 
As above, ^'j^ ^ define a linear isomorphism between H^ o and ^ using 

and the same i?jv,_D to define its converse explicitly. 



'AT 



3.2 Vector field on Bn^d lifted from the shape gradient 

Let VF be a shape, then for each a G Cj ([0, 1] , E) and for any Af e a ([0, 1]), 
VF associate to M an element VF (a) (M) g T^E. Consider now a e SJvd 
and ((Mo,o, • • • , Mo,d) , • • • , {Mn,o, Mn,d)) = Sn (a), then M^.q = Mn.d- 
We 

Tn,f (a) = ((VF (a) (M,,,) , j = . . . i?) , i = . . . i?) 
This representes the samphng of the deformation of the curve impUed by the 
shape gradient VF to a. It is not difficuh to sec that Tn.f (a) G T^jwia) ((^'^)^) • 

Theorem 1. To each shape gradient VF the map B^^jjoTn.f associate a vector 
field on Hjsi^D which correspond to a vector field Vp on B^j jj through jj. 

This theorem allows to interpret gradient descent method for shape optimiza- 
tion as a algorithm for integrating a vector field in a finite dimensional space. 
From this point of view, gradient descent method correspond to the most naive 
method to integrate this vector field, namely the Euler method. Clearly, this 
approach suggests to use better algorithm for vector field integration. 

3.3 Geometry of the vector field and local extrema of shape cost 
functional 

Proposition 9. Let a £ ([0, 1] ,E) be such that V F (a) = 0, i.e. VF (a) (M) = 

for all M £ a ([0,1]) and let 7 € Bm,d such that J {j^) — {j^) for 

z e {0, . . . , ND}, i.e.-i^ ^N,D {b-^]u ° (a)) , then VF (7) (7 (^)) = 

for i G {0, . . . , N and then Vp {^B~^^jj o Sn (a)^ = 0. It is to say that a local 
extremum of F induces a local extremum of its restriction to Bn.d o,nd that this 
extremum is "lifted" on a singularity of the vector field Vp on Hn,d- 

Proof. The deformation curve of 7 induces by the gradient of F vanishes at at 
least [D + 1) [N + 1) points, but it is a "Bezier curve" of degree D, so it a zero 
polynomial. So, its control polygon is reduce to the origin and then B^^j^oSn (ch) 
is a singularity of Vp- 

In fact, the vector field Vp is associated to the gradient of the function 
F o ^N,D- This is an heavy constrain on the vector field. For instance, it is easy 
easy to see that P is an attractive singularity of Vp if and only if 'Fn,d (P) is a 
local minimum of F | Hn d ■ 

3.4 Application to a problem of images segmentation 

In this section, we sketch an application to a problem of images segmentation. It 
is a problem of omnidirectional vision. Previous methods tried with some success 
but does not allowed a full real time treatment. There all based on snake-like 
algorithms (see [9]). The gradient use to detect edges is a classical one based 



on a Canny filter and is combined with a balloon force. The best previously 
known method is such that propagation of the contour were done using the fast 
marching algorithm for level set method. This a typical formulation of image 
segmentation as a shape optimization problem. In [S], we use piecewize Bezier 
curves to contour propagation and achieve a very fast segmentation algorithm 
allowing real time treatment even with sequential algorithm (no use of paral- 
lelism or special hardware architecture) on a embedded system. 




It is very interesting to see that, with few algorithmic modification, it is also 
possible to treat change of topology, i.e. curves with several connected compo- 
nents as it is shown in the following figure. 
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